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ABSTRACT 



A divergent heat budget equation which included the 
effects of surface heat flux, horizontal and vertical ad- 
vection, and horizontal divergence on the near-surface heat 
content was used to examine the role of thermal advection in 
the upper 250m of the water column at Ocean Weather Station 
NOVEMBER. This station is located on the southern boundary 
of the transition zone separating the Subarctic water mass 
from the Subtropic water mass. Values for horizontal thermal 
advection changes were computed over the period 1968-1970. 
Pulse-like periods of cool advection were associated with 
periods of reduced salinities suggesting these were intrusions 
of Subarctic water. Over the nine-year period of analysis, 
these intrusions had a periodicity of 7 to 8 months with a 
duration of 3 to 4.5 months. It is suggested these wave-like 
intrusions along the Subtropic front are the result of the 
passage of non-dispers ive baroclinic Rossby waves. 
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I. INTRODUCTION 



A. GENERAL 

The heat content of the upper layers of the ocean has 
recently been of interest, particularly with respect to the 
factors or forces which alter the heat state. Changes in 
the heat content result in changes in the temperature struc- 
ture and, with that, changes in the density distribution 
and oceanic layering. A knowledge of these factors and an 
attempt to predict and quantify their contributions is of 
importance to the United States Navy, especially with regard 
to anti-submarine warfare and its strong dependence on sound 
speed profiles which readily respond to alterations of the 
thermohaline structure. 

B. THE HEAT BALANCE 

An influential process in the alteration of the thermal 
structure begins at the air-sea interface with insolation, 
the rate at which short-wave energy from the sun enters the 
sea. This factor, combined with the effects of back radia- 
tion, evaporation, condensation, and conduction, produces a 
net surface heat exchange across the air-sea interface. In 
shallow water the effect of this net heat exchange is dis- 
proportionately large; in deep water its role as a significant 
process can vary with respect to other forces. 
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Redistribution of heat can occur via convective and 
mechanical mixing, the former primarily a seasonal process 
responding to the air-sea temperature differences, the latter 
responding to the effect of wind and waves. Both are 
significant in the layer above the main thermocline. 

Advection, vertical or horizontal, is a process of 
transport of an oceanic property (e.g., temperature) solely 
by the mass motion of the ocean. Its importance as a major 
heat contributor has been uncertain. Studies as recently as 
Thorne [1974] have computed horizontal thermal advection as 
the difference between the surface heat exchange (Q n ) and 
the local change in heat content ( 3 H/ 8 t). However, no cor- 
relations between the changes in horizontal advection and 
those of net surface heat flux, atmospheric pressure, or 
salinity were undertaken to support the validity of the model. 
It appears from this work and others [Tabata, 1961; Clark, 
1967; Bathen, 1970] that horizontal advection does not com- 
pletely account for the difference between local change in 
heat and net surface exchange; additional processes must be 
considered to account for this difference. 

Emery [1976] assumed that, in addition to Q n and hori- 
zontal advection (- 1 /^* V^H) , vertical advection and the 
compensatory action of horizontal divergence/convergence 
were major heat processes in the upper layer of the ocean. 
Correlation studies which included these additional terms 
suggested that they played a significant role in the heat 
budget of the upper ocean. Emery developed a divergent heat 
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budget (DHB) equation to model the correlation of these 
terms with the heat content of the upper layers. However, 
the contribution of each of the DHB terms to the heat balance 
was not specified nor were they correlated with changes in 
physical properties, e.g., salinity or oxygen. 

C. OBJECTIVES 

The purpose of this study was to complete and extend the 
findings of these earlier investigations. Specifically, the 
objectives of this research were: 

• to quantify each expression in the DHB equation 

• to determine the following trends 

- mean annual 

- yearly by month 

- anomalies 

• and in particular, to examine the significance of 
horizontal advection in changing the heat content of 
the upper ocean and to determine the validity of the 
results by correlation with other oceanic parameters. 

D. METHODOLOGY 

In order to measure .horizontal thermal advection 
(-t^h* VftH) , one must either possess a large horizontal data 
field of temperature and velocity to calculate the advective 
term directly, or deal with temperature and salinity measure- 
ments from a point source which necessitates accounting for 
all the expressions in the heat balance and attributing the 
residual to horizontal advection. In either case, it is 
important to validate the horizontal advective heat process 
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by comparing changes in the thermal structure with changes 
in a physical property such as salinity. Such a thermohaline 
re-structuring can be the result of water mass intrusions. 

Therefore, in order to realize the objectives of this 
research, a long-time data record is needed as well as a 
location that is readily influenced by water mass intrusions. 
Ocean Weather Station NOVEMBER (OWS N) affords such an oppor- 
tunity, since it offers both a long period of data and is 
strategically positioned at the transition zone between the 
Subarctic and Subtropical water masses. 

The Subarctic water mass is a cool, low-salinity body of 
water occupying the upper 300m of the Pacific Ocean, generally 
north of 40°N. Subtropical water is characterized by a 
shallow, warm, salty layer in the vicinity of 32°N. The 
transition zone between these water masses is characterized 
by sharp horizontal heat and salinity gradients [Roden, 1976] . 
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I I . THEORY 



A. GENERAL 

In his examination of the role of vertical motion in 
the heat budget of the upper 250m of the ocean, Emery [1976] 
found changes in the heat content (3H/3t) at OWS N could 
only be moderately explained by changes in heat across the 
air-sea interface (Q n ) . However, heat content changes could 
be explained far more accurately when two additional pro- 
cesses, horizontal divergence and vertical advection, were 
incorporated with the surface heat flux. To establish this 
correlation, he developed a divergent heat budget (DHB) 
equation which included these processes. 

The major assumptions in the development of the DHB 
equation are: 

• Molecular diffusion and radiative heat flux out of the 
column are neglected. 

• The field of motion is comprised of a steady mean flow 
upon which are superimposed internal Rossby waves with 
periods from months to years. 

• Temperature changes at the bottom of the water column 
can be attributed mainly to vertical advection of the 
mean thermal gradient at that depth. Since the base 
of the column is always beneath the main thermocline 
at OWS N, turbulent energy at this depth is small; 
therefore, temperature changes due to turbulent mixing 
can be considered negligible. 

Starting with the equation for the conservation of heat 
and by incorporating the above assumptions along with a 
Reynolds decomposition to the temperature and velocity 
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expressions, Emery arrived at the divergent heat budget 
equation, 



3H/3t - -jji- [H - P C p DT D ] + • V h H = Q n (1) 



1 



2 



3 



4 



where term 1 is the local change in heat content of a column 
of water of depth D; term 2 is the heat change due to hori- 



through the bottom of the column; term 3, the change due to 
horizontal advection; and term 4, the resultant gain or loss 
of heat across the air-sea interface. Term 2 may be combined 
as shown and expressed as H^ v , the net change in heat content 
due to near-surface divergence/convergence and the compensatory 
upwelling/downwelling mechanism. 

Q n is comprised of the following major terms, 



zontal divergence (W^H/D) and vertical advection 




Q n " Q s " Q b ' Q e Q h 



( 2 ) 



where 



Q s = the rate at which short-wave energy from the 
sun enters the sea. 



= the net long-wave back radiation from the sea 
surface . 



0 = the net latent heat transfer (includes 

x e 

evaporation and condensation) . 



0, = the net sensible heat transfer (conduction) . 

x h 
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H represents the heat content per unit surface area, 
evaluated for the upper 250m, and is computed by 



250 

H = 0 / pCpT dz 



(3) 



where p is the density and 
constant pressure. 

Wp, the mean vertical 



C is the 
P 



velocity , 



specific heat 



is computed by 



at 



W 



D 



~ 3T D/3t 

9T D/3z 



(4) 



which assumes that temperature changes at D are due mainly 
to vertical advection of the mean thermal gradient at D. 

^D/3t re P resents the monthly temperature change at 250m and 
3 Td/ 3 Z is the mean vertical temperature gradient at 250m. 

B. SIGN CONVENTION 

The following sign convention is explained to aid in 
later physical interpretation of the results (refer Figure 1) . 
Taking z positive downward and rearranging equation (1) to 
solve for the change in heat content in the water column, 



3H/3t = Q n + W d /D [H - pC p DT D ] - \ • V h H. (5) 

12 3 

1 . Net Surface Heat Exchange, Q n 

A positive value for Q n indicates that the column 
of water gains heat as the net result of the various heat 
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+ y (North) 



Surface 



► + x (East) 






250 m 



Figure 1: Direction of positive (heating) values 

for the DHB terms. 
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exchange processes acting at the air-ocean boundary 
[equation (2)]; negative values indicate that the ocean 
releases heat to the atmosphere via these processes. 

2. Divergent Heat Change , H^ iv = W^/D [ H - pC^DT^] 

a. [H - pCpDTp] is always positive since it 
represents pCp (T - T25Q ) where T is the average temperature 
throughout the column. 

b. From equation (4), can be either positive or 
negative. Examination of the records shows that at OWS N 
aTp/az is always negative at D = 250m. If 3T^/3t were 
positive, i.e., a heating process, then a positive results 
(downward vertical velocity) . A negative (upwelling) 
represents cooling. It follows then that term 2 of equation 
(5) responds directly to the sign of W^. Therefore, the 
cooling process is characterized by an upward vertical 
velocity (W^ negative) , accompanied by an out\^ard flow across 
the vertical sides of the column (divergence) in accordance 
with the law of continuity. Downward vertical velocity and 
inward flow are associated with heating. 

3. Horizontal Advection (-V^ * V.H) 

Expanding the horizontal advection term into its 
component parts yields 

-V h • V h H = -(ui + vj) • (3H/3x 1 + 3H/3y j). (6) 
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In the vicinity of OWS N isotherms are generally zonally 
distributed, hence horizontal advection reduces to(-v3H/3y). 
The north-south temperature gradient, 3T/3y is generally 
negative and can be appreciable when the southern boundary 
of the Subarctic/Subtropic transition zone migrates to the 
position of OWS N [Roden, 1976]. Positive horizontal advec- 
tion indicates that heat is carried into the column generally 
under conditions of a northward flowing current. A cooling 
process is generally in response to a flow from the north. 
When horizontal advection is due to the passage of an eddy, 
the flow direction cannot be readily inferred. Whether heat 
is gained or lost is dependent on the sign of the horizontal 
temperature gradient, i.e., whether a warm core or a cold core 
eddy is involved in the advection. 



% 
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III. ANALYSIS OF DATA 



A. SOURCE OF DATA 

Mechanical bathythermograph (MBT) data have been taken 
two to four times daily at OWS N since 1947. Monthly means 
of temperature at 5m intervals have been produced by Ballis 
[1973] from this data file. Significant gaps exist in these 
data for the earlier years, limiting a detailed analysis to 
the nine years, 1962-1970. In addition, MBT casts were 
limited to 135m prior to 1962, a depth too shallow for use 
in this study as the mixed layer is often slightly deeper 
than this depth, especially in winter. 

Within the 1962-1970 MBT time series, minor data gaps 
were eliminated by incorporating temperature data from the 
daily hydrocasts taken at OWS N since 1964 and on file at 
the National Oceanographic Data Center (NODC) . Although the 
monthly means calculated from the hydrocast data are based 
on fewer observations than those from the MBT data, no sig- 
nificant difference was detected between either data set. 
Standard deviations were similar and for those months when 
both sets of data existed, little difference in mean values 
was noted. 

A salinity time series was created from the NODC hydro- 
cast data. Data gaps limited this series to the period 1968- 
1970 since only a few short records existed prior to 1968. 
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Monthly means and anomalies of the net surface heat 
exchange (Q ) were taken from the report by Rabe and Bourke 
[1975] which examined the various surface heat exchange pro- 
cesses at OWS N. These data are based on three-hourly 
observations from whence mean daily values for. each month 
were calculated from standard bulk flux formulae. The 
dominant heat flux term, net solar insolation (Q ) , was 
calculated from the improved version of the Q s equation 
developed by Seckel [1973] . 

B. ANALYSIS TECHNIQUE 

Of the terms in the divergent heat exchange equation (1) , 
the primary focus of this report was directed towards the 
horizontal advection term (-V^ * V^H) . This is a residual 
term computed from the difference between the observed change 
in heat content (9H/3t) and the computed values of surface 
heat flux (Q n ) and divergent heat flux (H^ v ) . This type of 
analysis affords the possibility of error accumulation in the 
residual term. However, it will later be shown that, in spite 
of the possibility of error accumulation, the large anomalies 
in horizontal advection are significantly greater than the 
potential error contribution. 

1 . Determination of Local Heat Change 

Changes in the heat content of the water column are 
evaluated in this study on a seasonal scale by analyzing the 
temperature profiles from the MBT data of Ballis [1973] and 
the hydrocast data of Rabe and Bourke [1975] . These monthly 



21 



observations can be considered, in most cases, to be centered 
about the 15th of the month, since most monthly observations 
were taken over the entire 30-day period. Some monthly means 
were skewed to either side of the 15th in response to the 
biased nature of the data collection, e.g., temperature data 
for May 1968 were collected between the 12th and the 31st, 
shifting the effective mean date for this month to 21 May. 

It then follows that the average monthly heat content of the 
upper 250m, H, can be determined from these mean temperature 
profiles and assigned the same effective date. The mean heat 
content of the column is found from 

H = Z^pCpT dz (7) 

where "pC^ is the product of the mean water density and 

_ 7 _ 1 

specific heat, assumed constant at 0.978 g-cal cm °C , 
and T is the mean temperature profile over the specified time 
interval . 

The average monthly heat content was computed for the 
period 1962-1970; 3H/3t was subsequently determined for the 
interval between successive H's. The interval, 3t, ranged 
from 23 to 52 days. Values are computed in langleys per day 
(ly day' 1 ) . 

2 . Net Surface Heat Exchange 

The monthly mean daily values of Q , Q g , , and 
were combined to determine the monthly mean daily contribution 
of net surface heat exchange, Q n> 
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is the expression that 



3. Heat Divergence (H^ v ) 

As previously discussed, H, 
r ' ’ div 

includes heat exchange due to horizontal divergence and 
vertical advection, where W D = t [ 3 ^ D / 3 ' 1 • Similar 

to the procedure used to calculate 3H/3t, 3T D /3t was computed 
from the mean temperature at 250m, i.e., Tp = over the 

specified time interval, the interval between effective dates. 
The temperature gradient (ST^/Sz) at 250m was computed from 
the mean temperature profiles between 150 and 250m. This 
value ranged from -0.043 to -0.069 °C m" ^ . 

4 . Example of Computation Process 

The following is an example of the analysis technique 
applied to May and June 1968 (see Table I) . 

a. Effective Date 

As shown in the listing of the daily distribution 
of MBT observations for each month [Ballis, 1973], MBT casts 
for May 1968 were made daily between the 12th and the 31st; 
no observations were taken prior to the 12th. Therefore, the 
effective data for the mean monthly temp, T^j , is 21 May. 
Examination of the following month's record indicates at least 
one lowering per day was taken every day of the month (the 
most common situation). Hence the effective date, T j une > is 
15 June. The time interval (3t) between T^ and T June is 
established as 25 days. 
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Table I. Extract from mean monthly temperature profiles prepared 
by Ballis , 1973. 
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b. Heat Content 



The mean monthly heat content was computed by 
integrating the mean monthly MBT profile [Ballis, 1973] .to 
250m at 5m intervals. From equation (7), 



H May = 401 > 000 g m " ca -l c m 2 (langleys = ly) 



H June * 418 ’ 600 ^ 



It follows that the 3H/3t for this 25-day period 
is a heat gain of 17,600 ly or an average daily heat gain of 
700 ly day ^ . 

c. Vertical Velocity 
Recall that 




9T 2 5 q/ 9 1 
3 ^250^ 3Z 



(1) For May 



3T 



_250 = [11.49- 16. 10]°C = _ 0i0461 o c m 
9 3 250 - 150m 



-1 



(2) For June 



3T 



250 = [11-35 - 15, 901° C . . 0 . 0455 °C m‘ 



250 - 150m 



Therefore, the effective vertical temperature 
gradient during the 25-day interval separating the average 
profiles is -0.0458°C m ^ . 
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(3) Change in Temperature at 250m 



3 T 



250 _ [11.35 - 11 . 49] °C „ nrtr ^ 0 „ , -1 

3 = = -0.0056C day 

°t 25 days 



Therefore, 



T7T [-0.0056] n j -1 

W n = = -0.122 m day 

u [-0.0458] 



d. Average Net Surface Heat Exchange, Q 



n 



Average daily Q 



n 



May : -2 24 ly day 



June: 227 ly day 



-1 



-1 



The average net surface heat exchange for the period 21 May- 
15 June is: 



Q . (224 ] [10 days) + 227 [15 days] = 22S-8 ly day 
n 25 days 



-1 



e . H, . r 
div 



The FT^ term is dependent upon the mean heat 
content in the column between the effective dates (H) and the 
mean temperature at 250m over this same period CT ^ ^ q ) : 



W 

H J - •• - — c K - S DpT 250 ] 



div 



250m 



D = 250m 



From 4.b. and Table I, 



H = 418 > 6QQ ^ + 4Q1 . > . P - ° - 0 = 409,850 ly ; 



-=r _ 11.49 + 11.35 

A 250 ? 



11 . 42°C 
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f. Horizontal Advection ( • V^FT ) 

Solve the remaining terms in equation (1) 
where D = 250m : 



• ? h • 7 h H - W - 5 n = V [ H • pC p D t d 1 



W 

D r TT 



- - 225.8 - ['2^1] 



250 



[409,850 - (0.97 8) (2.5 x 10 4 cm) (11 . 42°C) ] 



538 ly day ^ (a heating process ) 



This 

heat column over 



is the average daily contribution to the 
the period from 21 May to 15 June 1968. 
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IV. RESULTS 



A. GENERAL 

Emery [1976] showed that the mean annual cycles of heat 
content within the mixed layer and within the upper 250m of 
the water column were significantly out-of-phase. This 
indicates that below the main thermocline, where turbulent 
kinetic energy is minimal, there are additional factors 
involved in the heat transfer process other than horizontal 
advection. Noting the fluctuations of the 14°C isotherm 
at both the short and long time scales, Emery assumed that 
vertical motion at 250m contributed significantly to the 
variability in the heat content and should be included in 
any conservation of heat equation. 

As illustrated in the preceding section, the magnitude 

of these vertical terms and the net surface heat exchange 

were computed from the monthly mean temperature profiles. 

The remaining contribution to the observed change in the heat 

content was attributed to horizontal advection, a term not 

computed specifically, but determined as a residual from the 

DHB equation. The mean local change in heat content (3H/9t) 

is plotted for the period 1962-1970 (Figure 2). A plot of 

the factors which contribute to the local heat change ( Q n , 

H.. , -V\ ♦ V, H ) are plotted for the same period in Figure 3. 

div’ h h J * 



28 





29 






Figure 3: Nine-year plot of the major 

contributors to the monthly heat 
change (3H/3t) at OWS N. 
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One immediately notices in Figure 2 the pulse-like occur- 
rences of large anomalies in the local heat change, especially 
the numerous events of large heat withdrawal, e.g., February 
1968 and April 1969. Comparing this time series with the 
corresponding values of the DHB terms of Figure 3, one notices 
that the pulse-like nature of the large anomalies in gH/9t, 
which is physically measured, cannot be satisfactorily accounted 
for by Q or H , . , since the former exhibits no pulse-like 
anomalies and those of the latter are small. It is horizontal 
advection, in fact, which must account for the large fluctua- 
tions in the gH/gt time series. Emery showed that upwelling, 
which is reflected in the H^ v term, operates to change the 
relationship between Q n and 9H/9t. However, this does not 
mean that H^ automatically acts as the major influence in 
the balance of heat. What Emery has done, in fact, in his 
examination of the role of vertical motion in the heat budget, 
is to make more reliable the approximation of the horizontal 
advection contribution. Therefore, of the DHB terms, the 
primary interest of this study lies in the role that horizontal 
thermal advection plays in altering the heat content of the 
upper layers of the water column from year to year. ^ Since 
-^h • V^H is a residually^determined expression which is 
potentially contaminated by measurement and modeling inaccu- 
racies, it is necessary to evaluate whether the magnitude and 
fluctuations of the term bear any resemblance to physical 
reality . 
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One possible approach is through water mass analysis, 
since the transport of heat is usually accompanied by the 
transport of other properties such as salinity or dissolved 
oxygen. Water mass analysis through temperature and salinity 
correlations is particularly useful at OWS N since it lies 
on the southern edge of the transition zone between the Sub- 
arctic and the Subtropical water masses, a frontal boundary 
distinguished by its sharp horizontal and vertical gradients. 
The position of OWS N is fortuitous since any significant in- 
trusion of cool, low-saline, Subarctic water should strongly 
influence the thermohaline structure of the upper 250m. 

Correlation of salinity changes with thermal advection 
changes at OWS N is necessarily limited to the three years, 
1968 through 1970, due to a lack of salinity data for the 
previous years. However, because of the good correlation 
between salinity and thermal advection during 1968-1970 
(shown later in this section) , one can surmise that thermal 
"events" of horizontal advection in the earlier years can be 
related to water mass intrusions. 

First, the data will be analyzed by examining the mean 
annual contribution from each of the terms in the DHB equa- 
tion, followed by a discussion of the major horizontal 
advective anomalies. 

B. MEAN ANNUAL CYCLE 

The mean annual cycle for each term in the DHB equation 
was calculated by averaging the mean monthly values from 



1962-1970 (Figures 4 to 7) . One must interpret the results 
with caution as the mean annual curves are based on a maximum 
of nine years of monthly averages, far too few to draw 
conclusive results. 

The mean annual cycle of 3H/3t (Figure 4) exhibits a 
broad general tendency for heating from April to September, 
followed by a cooling trend for the remainder of the annual 
cycle. The actual curve is a series of oscillations fluctu- 
ating about this trend by 200 to 300 ly day - '*'. However, of 
the component terms in the DHB equation, only Q n displays an 
annual cycle and is therefore the term which creates the 
cyclic trend in 3H/3t. The other terms appear to fluctuate 
as pulses and serve to enhance or to suppress the trend in 
surface heat flux. 

The pulse- like nature of the horizontal advection term 
appears non-seasonal and is noticeably larger in most cases 
than in its contribution to 3H/3t. If, as anticipated, 

its oscillatory nature is due to a water mass intrusion, 
the thermohaline structure of the environment should be 
significantly altered. Hence, salinity changes during any 
large advective cycle should also exhibit a pulse-like nature. 
This hypothesis is supported by Hansen [1973] who observed 
that surface salinity changes at OWS N were pulse- like occur- 
rences displaying neither a seasonal trend nor a correlation 
with any seasonal process, e.g., evaporation minus precipita- 
tion . 
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LANGLEYS 




Figure 4: Mean annual trend in the monthly changes 

in heat content (9H/3t) at OWS N. Bars 
indicate one standard deviation (±lo) 
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LANGLEYS 
PER DAY 




Figure 5: Mean annual trend in the horizontal 

advection contribution (-V^ • V^H) 
to local heat change (3H/9t) at 
OWS N. Bars indicate one standard 
deviation (±lo) . 
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LANGLEYS 
PER DAY 




Figure 6: Mean annual trend in horizontal 

divergence plus vertical advection 
contribution (H , , ) to the local 
heat change (3H/3t) at OWS N. 

Bars indicate one standard deviation 
(±lo) . 
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LANGLEYS 
PER DAY 




Figure 7: Mean annual trend of net surface heat 

exchange contribution (Q n ) to local 
heat change (3H/3t) at OWS N. Bars 
indicate one standard deviation (±la) . 
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C. HORIZONTAL ADVECTION ANOMALIES 



The years 1968-1970 were highlighted by five distinct 
anomalies or events in horizontal advection. A detailed 
examination of each event is presented in Appendix A. 

The initial phase of each of these irregularities was 
a sharp cooling trend followed by a substantial recovery 
period of warm advection. Rapid cooling of the water column 
occurred over a period of 25 to 30 days causing the heat 
reduction by horizontal advection to increase from near zero 
values to approximately 300 to 850 ly day ‘ L at the peak of 
the cold water intrusion. Thus, the daily change of heat 
loss during these cooling periods was approximately 10-23 ly 
day ^ . 

A rebound to initial conditions occurred following these 
cooling peaks over approximately the same time period as the 
cooling trend. Typically, the return to mean conditions was 
followed by an overshoot of heating, peaking at 300 to 
500 ly day ^ , eventually returning to mean conditions about 
3 to 4.5 months after the intrusion of the low- temperature 
water mass. 

Salinity changes during each of these events were signif- 
icant at the surface and throughout the 250m water column 
(see Table II, Appendix A). Each thermal advection event 
was accompanied by a concurrent salinity change, i.e., the 
initial cooling trend was accompanied by a large salinity 
decrease throughout the column, while the rebound phase of 
warm advection was characterized by salinity increases. 
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The average changes at the surface, 50m, 150m, and 200m, 
during cooling trends were -0.21, -0.18, -0.17, -0.15°/ oo , 
respectively. During the rebound phase mean salinity increases 
were 0.23, 0.16, 0.13, 0.10°/ O o* 

Since each horizontal advective cooling trend was char- 
acterized by a concurrent and significant decrease in the 
mean salinity and each advective heating trend was accompanied 
by an equally sharp salinity increase, it can be concluded 
that the major cause for these thermohaline fluctuations was 
the intrusion of the Subarctic and the Subtropic water masses. 

Also note that in the nine-year time series (Figure 3) 
horizontal advection (-V^ ♦ V^H) and H^ v do not necessarily 
fluctuate in phase. For example, refer to Events 1 to 5 in 
Figure 8. Note that in Events 1, 2, 3, 8 5 the H^ v term 
follows the trend of horizontal advection; in Event 4 H^ v 
is oppositely directed. 

An interesting phenomenon was noted in the salinity 

changes at 200m during these events. Events 1, 2, 3, 8 5 

are periods of cool horizontal thermal advection, i.e., 

intrusions of low-salinity Subarctic water. Simultaneously 

with these Subarctic intrusions, H , . indicates the occurrence 

div 

of upwelling. During these trends, the salinity decrease at 
200m averages -0. 12°/ 0 o - However, in Event 4, where the 
term indicates that downwelling was occurring, the salinity 
change at 200m is a very large -0.34 o / oo . A possible explan- 
ation for these salinity differences is that during each of 
these intrusions of Subarctic water, the temperature and 
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Figure 8: Major events in horizontal advection (-V7 • V, H) for 1968-1970. 



salinity in the upper layers are lowered. During coincident 
upwelling periods (H diy negative) , the salinity decrease due 
to advection is suppressed by the introduction of higher 
salinity water from the lower portion of the Subarctic halo- 
cline (Figure 9). However, during Event 4 the downwelling 
action enhances the salinity reduction due to advection by 
mixing downward the lower salinity values of the upper 
halocline . 

D. BAROCLINIC ROSSBY WAVES 

These thermohaline changes suggest the intrusion and 
subsequent recession of the Subarctic water mass at OWS N in 
a pulse-like, non-seasonal manner. This translation of the 
boundary may be explained only partially by examination of 
surface wind stress since the resultant Ekman dynamics cannot 
fully account for the temperature and salinity fluctuations 
between the main thermocline and 250m. 

The non-seasonal nature of the horizontal advection 
anomalies indicates that perhaps the driving mechanism itself 
is non-seasonal. An examination of atmospheric data from 
1968-1970 reveals no unusual atmospheric pressure anomalies 
or significant weather deviations during the time of the large 
advective anomalies. What is apparent, however, about these 
pulse-like cooling/heating trends is their 7- to 8-month 
periodicity -- an indication perhaps that the southern boundary 
of the Subarctic/Subtropic transition zone might be responding 
to the passage of non-dispers ive baroclinic Rossby waves, as 
suggested by Bernstein and White [1974]. 
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Figure 9: Typical temperature, salinity, and density structure in 

the Subarctic Pacific Ocean [after Tully, 1963]. 



In their study of the time and length scales of baroclinic 
eddies in the Central North Pacific Ocean, Bernstein and White 
examined the subsurface temperature fluctuations at OWS N 
using the MBT data compiled by Ballis [1973] . They determined 
that the large-amplitude fluctuations of the average monthly 
temperature at 200m, which was below any seasonal influence, 
could not be due to effects of internal waves or variations 
in the geostrophic flow. They noted a long-period fluctuation 
of six years upon which smaller variations occurred with 
periods of two to six months. They attributed these scales 
of motion to baroclinic, non-dispersive, Rossby waves. How- 
ever, they were unable to localize the dominant period of 
maximum energy in the spectral density functions for the sub- 
surface temperature profiles as their data record length (the 
same as used in this report) limited the maximum resolution 
to about six months. Much energy is undoubtedly contained in 
longer period waves. From his study of spectra from tempera- 
ture records from the ocean weather ships in the Atlantic, 

Gill [1975] showed that the spectra are different at each 
station, but that large energy levels existed at periods 
greater than six months, especially for stations near frontal 
boundaries. It is suggested, then, that the 7- to 8-month 
occurrences of large advective cooling/heating, evident in 
this analysis, are the results of baroclinic Rossby waves. 

E . SUMMARY 

It is apparent that the effect of these wave-like eddies 
on the temperature and salinity structure near water mass 



43 



boundaries can be substantial. Each of the large advective 
anomalies was characterized by an initial cooling trend, an 
indication that the southern edge of the transition zone at 
OWS N was being pushed southward. The thermal and salinity 
structure of the water column reacted accordingly to the 
intrusion of fresh, cool water from the north. As the eddy 
progressed westward, it allowed the transition zone boundary 
to "rebound" towards the north wherein the thermal and salin- 
ity structure responded to the return of the warm, saline 
subtropical water mass. 
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V. CONCLUSIONS 



This study investigated the effects of horizontal thermal 
advection in the upper 250m at OWS N from 1962-1970. 

In accordance with the divergent heat budget (DHB) equa- 
tion, heat changes that were not due to the combined effects 
of net surface exchange, vertical advection, and horizontal 
divergence were attributed to horizontal thermal advection 
in maintaining the balance of heat in a column of water. 

The investigation concentrated on the period 1968-1970 
and centered on water mass analysis and its relation to the 
large anomalies in horizontal advection. The results of 
this research show: 

• The divergent heat budget equation is a valid approach 
to the determination of horizontal advection. 

• Q n , the net surface heat exchange, exhibits an annual 
cycle and is therefore the term which creates the 
cyclic trend in aH/at. The other DHB expressions 
fluctuate as pulses and serve to enhance or suppress 
the trend in surface heat flux. 

• The large pulses in 3H/at are the result of large 
pulses of horizontal thermal advection and, to a 
small degree, are influenced by the pulses of 

• During the years 1968-1970, there were five significant 
horizontal advective anomalies, each distinguished by 
an initial cooling trend in conjunction with a sharp 
salinity decrease, followed by an advective warming 
period and salinity increase throughout the entire 
250m column. These large advective cycles coincided 
with the largest of the changes in mean monthly salinity. 

• The change in salinity is in response to translation 
of the Subtropic/Subarctic water mass boundary zone. 

Cold horizontal advection is a result of the intrusion 
of Subarctic water while advective warming is indicative 
of a return of the Subtropical water mass. 
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• Horizontal advection and H^iy do not necessarily 
fluctuate in phase. A period of upwelling (H(j lV 
negative) serves to suppress the salinity change at 
200m due to horizontal advection by mixing upward 
the higher salinities from the lower portion of the 
intruding Subarctic halocline. Conversely, down- 
welling (H^y positive) enhances the salinity decrease 
due to horizontal advection. 

• The pulse-like nature of these five events is suggested 
to be in response to a non-dispers ive baroclinic Rossby 
disturbance whose wave-like characteristics are respon- 
sible for the depression of the southern edge of the 
Subtropic/Subarctic transition zone into OWS N. 
Intrusions of Subarctic water occur approximately every 
7 to 8 months. The duration of each intrusion is about 
3 to 4.5 months . 

• Other horizontal advective anomalies occurred in the 
period 1962-1967, but remain unsubstantiated by any 
changes in oceanic physical properties. In particular, 
a distinct lack of documented salinity data for this 
period prevented the examination of any thermohaline 
response to horizontal advection. However, these 
anomalies are, as in 1968-1970, non-seasonal in occur- 
rence and exhibit the same characteristic period and 
duration. 
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APPENDIX A 



HORIZONTAL ADVECTION EVENTS 

The details of each horizontal advection event are pre- 
sented in Figure 10. Each event is divided into a cooling 
trend (A) and a heating trend (B) . Each of these trends is 
associated with a salinity change throughout the column, 
values for which are incremented at the surface, 50m, 150m, 
and 200m and are listed in Table II. Note that cool advec 
tion is always associated with a reduction in salinity and 
that salinity increases during periods of advective warming 
The large change in salinity at 200m during Event 4 occurs 

because H , . is out of phase with horizontal advection. 
div 
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l’igure 10: Major events in horizontal advection • V, II) for 1968-1970. 



Table II: Salinity changes (°/oo) during 

horizontal advection anomalies 







Sal inity 


Change 






Event 




Surface 


50m 


150m 


200m 


1A 




-0.26 


-0.19 


-0.20 


-0.18 


2A 




-0.14 


-0.18 


-0.11 


-0.15 


3A 


-Vv 


-0.11 


-0.13 


-0.20 


-0.10 




Cooling 










4A • 




-0.14 


-0.15 


-0.10 


-0.34 


5A 




-0.28 


-0.17 


-0.12 


-0.06 


IB 




+ 0.25 


+ 0.19 


+ 0.14 


+ 0 .23 


2B 




0.24 


0.13 


0.05 


0.01 


3B 


-Vv 1 


0.26 


0.20 


0.30 


0.22 




heating 








0.08 


4B 




0.30 


0.27 


0.11 


5B 




0.10 


0.12 


0.10 


0.05 
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